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f.: Abstract 

;_^ ' We investigate the effects of newly measured /3-decay half-hves on r-process nucleosynthesis. 
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C^ ' These new rates were determined by recent experiments at the radioactive isotope beam factory 



facility in the RIKEN Nishina Center. We adopt an r-process nucleosynthesis environment based 
on a magnetohydrodynamic supernova explosion model that includes strong magnetic fields and 
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^vq ' rapid rotation of the progenitor. A number of the new /3-decay rates are for nuclei on or near the r- 

in 

"^ • process path, and hence, they affect the nucleosynthesis yields and timescale of the r-process. The 

^^ i main effect of the newly measured /3-decay half-lives is an enhancement in the calculated abundance 

of isotopes with mass number A = 110 - 120 relative to calculated abundances based upon /3-decay 

k>i ■ rates estimated with the finite-range droplet mass model. This effect slightly alleviates, but does 

C^ . not fully explain, the tendency of r-process models to underproduce isotopes with A = 110 - 120 

compared to the solar-system r-process abundances. 
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Introduction. Rapid neutron-capture (r)-process nucleosynthesis is responsible for the 
origin of approximately half of the elements heavier than iron and is the only means of 
producing the naturally occurring radioactive heavy actinide elements such as Th and U. In 
spite of more than a half century of study and observational progress [e.g., |l|, however, the 
astrophysical sites for r-process nucleosynthesis have still not been unambiguously identified 
(for recent reviews see, e.g., [21, yf). Although many candidate sites have been proposed and 
supernovae appear to be well suited as the r-process site 4!], up till now there has been no 
consensus as to the correct astrophysical model. 



e_ astronomical environment, the phys- 
It is evident that the r-process oc- 



Notwithstanding the difficulties in finding a suitab 
ical conditions for the r-process are well constrained 
curs via a sequence of near equilibrium rapid neutron captures and photo-neutron emission 
reactions far on the neutron-rich side of stability. This equilibrium is established with a 
maximum abundance strongly peaked on one or two isotopes far from stability. The relative 
abundance of r-process elements is then determined by the relative /3-decay rates along this 
r-process path., i.e., slower /3-decay lifetimes result in higher abundances. At least part of 
the reason for the difficulty in finding the astrophysical site for the r-process stems from the 
fact that it lies so far from the region of stable isotopes where there is little experimental 
data on /3-decay rates. 

In this context, it is of particular interest that /3-decay half-lives of 38 neutron-rich iso- 
topes including looKr, lo^-iosg^^ loe-iosy^ 108-1102^, m.nsNb, "^-hsmo, and 116,117^^ ^^^^^ 

been measured J6| at the recently commissioned radioactive isotope beam factory (RIBF) 
facility at the RIKEN Nishina Center. A secondary beam, containing the neutron-rich nuclei 
of interest, was produced by infiight fission of a 345 MeV/nucleon ^ss-q i^g^j^i in a Be target. 
The nuclei in the secondary beam were identified on an event-by-event basis. Their atomic 
numbers were determined from the energy loss in an ionization chamber, and the charge to 
mass ratio was determined by combining the projectile time-of-fiight and magnetic rigidity 
measurements. The nuclei were then implanted in a silicon-strip detector where the /3-decay 
lifetimes could be measured. Most of the measured lifetimes are an improvement on existing 

n 

measurements and a number of them [6| were measured for the first time. 

Many of these isotopes are near or directly on the r-process path. As such, they are 
of particular interest as they determine the /3-fiow toward the important r-process peak 
ai A = 130. Thus, they regulate the ability of models for the r-process to form heavier 



elements [7|. It is of particular interest, therefore, to examine the impact of these new rates 
on r-process models. 

The magnetohydrodynamic (MHD) supernova explosion model is of particular interest in 
the present study as it provides a good fit to the observed r-process abundances and avoids 
some of the problems |7|, and references therein] associated with other paradigms such as 
the neutrino driven supernova wind scenario |8|. Robust r-process nucleosynthesis appears 
to be a natural consequence of MHD mechanisms for core-collapse supernova explosions 
involving strong magnetic fields and rapid rotation. These models are also consistent with 
the observed jets and asymmetry of core-collapse supernovae. The detailed hydrodynamic 
properties of the MHD supernova model and the input nuclear physics are described below. 

Magnetohydrodynamic Models. The main features of the MHD supernova model employed 
here are described in Ref. |9i] and need not be repeated here. However, we emphasize that 
although this calculation is performed in the context of a particular r-process paradigm, this 
model is selected because it matches the required conditions of the r-process environment 
(timescale, neutron density, temperature, entropy, electron fraction, etc.). Thus, the results 
presented here should be similar to those obtained from more generic r-process models [e.g., 
7|. 
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We utilize the jet-like explosion (model 4 in [9[) based upon a two-dimensional mag- 
netohydrodynamic simulation. The ejecta were evolved with 23 tracers to describe the 
evolution of thermodynamic state variables. These were then post-processed to obtain the 
nucleosynthesis yields. Figure [1] shows the evolution of entropy, measured in units of the 
Boltzmann constant ks, and electron fraction, Fg, after nuclear statistical equilibrium (NSE) 
is achieved for three representative tracer particles together with temperature and density. 
Hydrodynamic quantities (e.g., entropy, temperature, and density, etc.) are obtained from 
the magnetohydrodynamic simulation, whereas the electron fractions are deduced from the 
nucleosynthesis calculation. We switch from the NSE abundances to a nucleosynthesis net- 
work calculation at T = 9 GK (~ 0.78 MeV). 

Nuclear Reaction Network. The nuclear reaction network utilized for the nucleosynthesis 
simulations has been described in detail elsewhere J9|, Il0| and need not be described in more 
detail here. The network consists of more than 4000 isotopes, including neutrons, protons. 



and heavy isotopes with atomic number Z < 100 (for detail, see Table 1 in Ref. 
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We also 



consider possible reactions related to the r-process involving two- and three-body reactions 
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FIG. 1. Thermodynamic properties for three representative tracer particles are shown as A {red 
solid lines), B (green dotted lines), and C {blue dashed lines). From top to bottom, each figure 
shows the time evolution of temperature T in GK, density p in g cm^'^, entropy per baryon S/k^ 
in Boltzmann constant k-Q, and electron fraction Yg, respectively. The evolution of each variable is 
plotted for 0.5 s from the time of the switch from NSE to a network calculation at T = 9 GK. 



or decay channels, and we include electron capture as well as positron capture and screening 
effects for all of the relevant charged particle reactions. Experimentally determined masses 



ll| and reaction rates are adopted if available. Otherwise, the theoretical predictions for 



nuclear masses, reaction rates, and /3-decays are obtained from the finite-range droplet model 
(FRDM) pj. 

In the present study, we perform r-process calculations based upon three different nu- 
clear reaction networks as summarized in Table [B These are extensions of the basic network 
described above. One of the networks utilizes only the FRDM theoretical rates from the 
REACLIB compilation [l3|. The other two (RIBF and RIBF-f) utilize the new experimental 
/3-decay half-lives of 38 neutron-rich isotopes from Kr to Tc and two versions of the theo- 
retical FRDM /3-decay rates for the other isotopes. The RIBF network replaces the FRDM 
rates with the new measured ones where possible. The network (RIBF+) is based on the 
RIBF network and FRDM rates with modified Q values for (n, 7) and reverse reactions given 

by: 

f 

Q - 0.3 [MeV] ( 97 < A < 103) 

Q + 0.5 [MeV] (104 <A< 107) (1) 

Q + 1.0 [MeV] (108 <A< 115) 

where Q is the theoretical Q value (in MeV) obtained from the FRDM [12] . The motivation 
for these modified Q values is that they lead to a better fit to the observed /?-decay lifetimes 
away from stability. 

Results. The nucleosynthesis calculations have been performed for 23 tracer particles 
from the MHD supernova model (e.g. Figure [1]) with the three different reaction networks 
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TABLE I. Reaction rates used in nuclear reaction networks 



/3-decay half-lives 



(n,7), (7, n) reaction 



FRDM 

RIBF 

RIBF+ 



REACLIBf. 
REACLIB + RIKEN^ 
REACLIB + RIKEN 



^ REACLIB : the REACLIB compilation fn 
^ RIKEN : experimental data given in Ref. 
^ Q+ modified : defined in Eq. ^. 



Ref. 
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REACLIB 
REACLIB 
REACLIB + Q+ modifiecP 



listed in Table HI In particular, the underproduction of isotopes near A = 120 becomes 
slightly less pronounced relative to predictions based upon the FRDM rates when the new 
measured rates are employed. Figure [2] illustrates the final abundances obtained from the 
three representative trajectories (A, B, and C) of Figure [T] and the three reaction networks 
employed. Note that as the A ~ 130 r-process peak is approached, the new rates begin to 
make a difference in isotopes with A ^ 90 - 120. 

Figure [3] shows the final integrated abundance distribution from all of the trajectories. 
These are compared to the solar system r-process abundance distribution [ij]. Here the 
effect of new rates becomes more apparent. There has been a recurrent conundrum in the 
r-process models in that they tend to underproduce nuclei with A ~ 120 [8]. One hope 
has been that the newly measured /3-decay rates in this mass region might shift the /3-flow 
equilibrium thereby filling in the low abundances near A ~ 120. Here, however, we see on 
Figure [3] that the abundances in the A = 110 - 120 region are only slightly enhanced. Thus, 
although the new rates provide a little assistance in enhancing the abundances near the 
valley, they do not alleviate this problem. 

This suggests that a further modification of the r-process paradigm is required. One 
possibility is suggested by this study. We note that the r-process fiow moved farther away 
from stability and proceeded faster in the RIBF+ network in which the (n, 7) Q values 
were systematically enhanced for isotopes with A = 104 - 115. This helped to fill in the 
abundances in the higher mass region with A > 115 (see Figure |3]). Thus, it is important 
to measure the masses and/or neutron separation energies in this region. For example, if 
the strength of the nuclear closed shell near the neutron magic number 82, and A = 120 

- 140 was systematically diminished, this would prevent the r-process path from bypassing 
the A ~ 120 nuclei as the (r2,,7) equilibrium shifts toward the more strongly bound nuclei 
along A^ = 82 the neutron closed shell. Whatever the explanation for the filling of ^4 ~ 110 

- 120, however, it is evident that the new /3-decay rates have provided new insight into the 
nucleosynthesis of the heavy nuclei in the r-process. 

At the very least, these results confirm that the final abundances in the mass region 
oi A ~ 110 — 120 are quite sensitive to the /?-decay half-lives of isotopes along the r- 
process path. Indeed, the astrophysics of the production of nuclei in the mass region of 
A ~ 110 — 120 is currently of considerable interest. It is presently thought that the lighter 
heavy elements with A < 120 observed in some ultra-metal-poor stars, may be the result of 
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FIG. 2. Upper, middle, and lower panels show final abundance distributions for the illustrative 
tracer particles A, B, and C from Figured) Each panel shows three different lines corresponding 
to the different reaction networks, i.e., FRDM (solid red lines), RIBF (dashed green lines), and 
RIBF+ (dotted blue lines). In the case of tracer A, the three lines are indistinguishable because 
the abundance distribution does not reach mass numbers A ~ 100 - 130 where differences in the 
networks appear. 



< 

> 

o 

D) 
O 



-5 -5 



-6 7 



— 1 1 1 1 1 1 r 



—I I I i_ 



FROM 

RIBF 

RIBF+ 




_l I I L 



- -6 



50 100 150 200 90 100 110 120 
mass number, A mass number, A 



FIG. 3. Integrated mass averaged total final abundance distributions of r-process elements from 
the adopted MHD supernova model (jet model 4 in Ref. [Qi]). Red solid, green dotted, and blue 
dashed lines correspond to results from using the FRDM (standard), RIBF, and RIBF+ rates, 
respectively. Abundances of solar system r-elements [ij] are represented by black dots with error 
bars. 



a new light-element primary process (LEPP) [15]. Our results suggest that nucleosynthesis 
in the LEPP should also be sensitive to the nuclear physics uncertainty from /3-decay rates 
in this region. Hence, further studies on both the nuclear physics and astrophysics of the 
synthesis of elements with A ~ 110 — 120 are warranted. 
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